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Abstract 

We present results for the structure function Fl for a gluon target having nonzero 
transverse momentum square at order a s . The results of double convolution (with 
respect to Bjorken variable x and the transverse momentum) of the perturbative 
part and the unintegrated gluon densities are compared with recent experimental 
data for Fl at low x values and with predictions of other approaches. 
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1 Introduction 



The basic information on the internal structure of nucleons is extracted from the process 
of deep inelastic (lepton-hadron) scattering (DIS). Its differential cross-section has the form: 



where F 2 (x, Q 2 ) and F L (x, Q 2 ) are the transverse and longitudinal structure functions (SF), 
respectively, and are the photon and the hadron 4-momentums and x = Q 2 / (2pq) with 

Q2 = _ q 2 > o 

The longitudinal SF Fl(x, Q 2 ) is a very sensitive QCD characteristic because it is equal to 
zero in the parton model with spin— 1/2 partons. Unfortunately, essentially at small values 
of x, the experimental extraction of Fl data is required a rather cumbersome procedure (see 
[l], ||, for example). Moreover, the perturbative QCD leads to some controversial results in 
the case of SF Fl- The next-to-leading order (NLO) corrections to the longitudinal coefficient 
function, which are large and negative at small x ||, |4j], need a resummation procedure [] 
that leads to coupling constant scale higher essentially then Q 2 (see || [|, [7J) Q 

Recently there have been important new data [pl|-[[18] of the longitudinal SF Fl,, which 



have probed the small- x region down to x ~ 1(T 2 . Moreover, the SF F L can be related at 
small x with SF F 2 and the derivation dF 2 /dln(Q 2 ) (see flI5|-|21|). In this way most precise 



predictions based on data of F 2 and dF 2 /d\n(Q ) (see [|L5[] and references therein) can be 



obtained for Fjj. These predictions can be considered as indirect 'experimental data' for Fl. 
In this paper for analysis of the above data we will use so called /^-factorization approach 



22 , |23| . 24] based on BFKL dynamics [12| (see also recent review |TT| and references therein) 



In the framework of the /^-factorization approach, a study of the longitudinal SF Fl has been 
done firstly in Ref. [25], where small x asymptotics of Fl has been evaluated using the BFKL 
results for the Mellin transform of the unintegrated gluon distribution and the longitudinal 
Wilson coefficient functions has been calculated analytically for the full perturbative series at 
asymptotically small x values. Since we want to analyze Fl data in a broader range at small 
x, we will use parameterizations of the unintegrated gluon distribution function & g (x, k'jj 
(see Section 3). 

The unintegrated gluon distribution & g (x, k\) (f g is the (integrated) gluon distribution 
in proton multiplied by x and kj_ is the transverse part of gluon 4-momentum k^) 

rQ 2 

f g (x,Q 2 ) = / dkl$ g (x,kl) (hereafter k 2 = -k]_) (1) 

is the basic dynamical quantity in the A>r-factorization approach It satisfies the BFKL 
equation 



1 Without a resummation the NLO approximation of SF Fl can be negative at low x and quite low Q 2 
values (see ||, ||). 

2 Note that at low x a similar property has been observed also in the approaches ||, ^, (see recent 
review [jllj and discussions therein), which based on Balitsky-Fadin-Kuraev-Lipatov (BFKL) dynamics jj^ j, 
where the leading ln(l/x) contributions are summed. 

3 In our previous analysis pBfl we shown that the property fc 2 = — k\ leads to the equality of the Bjorken 
x value in the standard renormalization-group approach and in the Sudakov one. 
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Notice that the integral is divergent at lower limit (at least, for some parameterizations 
of $ g (x, kj_)) and so it leads to the necessity to consider the difference f g (x, Q 2 ) — f g (x, Qq) 
with some nonzero Qq (see discussions in |26| 



i.e. 



f g (x,Q 2 ) = f 9 (x,Ql)+ f Q o dk 2 ± $ g (x,k] 



(2) 



Then, in the /^-factorization the SF F 2 ^(x, Q 2 ) are driven at small x primarily by gluons 
and are related in the following way to the unintegrated distribution § g (x, k'jj: 



F 2 , L (x,Q 2 ) = f'-fdkl £ e 2 .Cl L {x/z,Q 2 ,mlkl)^ g {z,kl), 

Jx Z J , 



(3) 



i=u,d,s,c 



where ef are charge squares of active quarks. 

The functions C 2j l(x, Q 2 > m i > can be regarded as SF of the off-shell gluons with 
virtuality k\ (hereafter we call them as hard structure functions Q). They are described by 
the sum of the quark boxes (and crossed boxes) diagram contribution to the photon-gluon 
interaction (see Fig. 1). 

The purpose of the paper is to give predictions for the longitudinal SF Fl(x,Q 2 ) based 
on the calculations of the hard SF C 2 L (x, Q 2 , m 2 , k\), given in our previous study [^6 



and several parameterizations of unintegrated gluon distributions (see JTTJ and references 
therein) . 

It is instructive to note that the diagrams shown in Fig. 1. are similar to those of the 
photon-photon scattering process. The corresponding QED contributions have been calcu- 
lated many years ago in Ref. p7| (see also the beautiful review in Ref. |28j). Our results 
have been calculated independently and they are in full agreement with Ref. |27| . Moreover, 
our results are in agreement with the corresponding integral representations for C 2L , given 
22, 25 1 and numerically with results of p9 |. However, we hope that our formulas which 



in 



are given in a simpler form could be useful for others. This simpler form for the hard SF C| L 
comes from using the relation between results based on non-sense, transverse and longitudinal 
gluon polarizations (see Eq. (|I3|) below) observed in for gauge-invariant sets of diagrams. 



The structure of this paper is as follows: in Section 2 we present the basic formulae of 
our approach. Section 3 contains the relations between SF Fl and F 2 and the derivative 
dF 2 /d In Q 2 , obtained in |19|, ^(| ^lj in the framework of Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) approach |30[ (i.e. in the collinear approximation: kj_ = 0). In Section 4 
we give the predictions for the structure function F^ for three cases of unintegrated gluon 
distributions. 



2 Basic formulae 



To begin with, we review shortly results of Ref. |26[ needed below in our investigations. 
The hadron part of the DIS spin-average lepton-hadron cross section can be represented 
in the form 

(4) 



in/ 



e^ u (q) F L (x,Q 2 ) +d IMU (q,p) F 2 (x,Q 2 ) 



4 by analogy with similar relations between cross-sections and hard cross-sections. 
5 Hereafter we consider only one-photon exchange approximation. 
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where 



and d^q.p) 



^ + 2a;fj h Vp.V 



4x| 

i/ 9 



2.1 Fey nman- gauge gluon polarization 



As it has been shown in Ref . |26| , it is very convenient to consider, as the first approximation, 
gluons having polarization tensor (hereafter the indices a and (3 are connected with gluons 
and fi and v are connected with photons)^: 



pa/3 _ _g a P 



(5) 



The tensor corresponds is equal to the standard choice of polarization matrix in the frame- 
work of collinear approximation. In a sense the case of polarization is equal to the standard 
DIS suggestions about parton properties, excepting their off-shell property. The polarization 
(|) gives the main contribution to the polarization tensor we are interested in (see below) 



BFKL 



ha h.0 



(6) 



which comes from the high energy (or kj) factorization prescription |22, 24] []. 

Contracting the photon projectors (connected with photon indices of diagrams on Fig.l.) 



p(i) 



\ 9 „ and P^=Az 2 ^ 



(here z = Q 2 / (2pq) is the corresponding Bjorken variable on parton level) with the parton 



tensor F*: 



= e, v {q) F p L (z,Q 2 ) + d, u (q,k) F?(z,Q 2 ), 



(7) 



we obtain at the parton level (i.e. for off-shell gluons having momentum fc M ), when C% l(z) 
F 2 P L (z,Q 2 )$ 



P 2 ■ Ci(x) 
P 2 ■ Cl{x) 



/(!) + 



2(3 2 



(2) 



Abx 2 f (1) + 



f 



;i + 2bx 2 ) 



f 



(2) 



jc . /( 2 ) + 4bx z p 2 ■ C g 2 



2 52 



(8) 
(9) 



6 In principle, we can use here more general cases of polarization tensor (for example, one is based on 
Landau or unitary gauge). The difference between them and (0) is ~ k a and/or ~ and, hence, it leads 
to zero contributions because Feynman diagrams on Fig.l are gauge invariant. 

7 We would like to note that the BFKL- like polarization tensor (||) is a particular case of so-called nonsense 
polarization of the particles in i-channel. The nonsense polarization makes the main contributions for cross 
sections in s-channel at s — ► oo (see, for example, Ref. jH] and references therein). The limit s — > oo 
corresponds to the small values of Bjorken variable x, that is just the range of our study. 

8 The hard SF CS L do not depend on the type of target, so we can replace z — ► x below. 
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where the normalization factor K, = a s (Q 2 ) ■ x, 

Pj$F ia , = lC-f {fi ,i = 1,2 

and a s {Q 2 ) = a s (Q 2 )/(47r), /f = 1 - Abx 2 , b = -k 2 /Q 2 = k 2 JQ 2 > 0, a = m 2 /Q 2 . 

Applying the projectors P^J to Feynman diagrams displayed in Fig.l, we obtain the 
following results 



./ 



(i) 



1 - 1 - 2x(l + b - 2a) ■ [1 - x(l + b + 2a)] ) ■ f\ 



-2(3 



+ {2a-b){\ -2a)x 2 ■ f 2 



f 



(2) 



(1 - (1 + b)x) - 2x[bx(l - (1 + b)x)(l + b - 2a) + a(3 2 ■ f x 



lx-f3 



+ bx 2 {\ - (1 + b)x){2a - b) -h 



(10) 



11) 



where 



(3 2 = l- 



4ax 



1 - (l + b)x) 



(12) 



and 



, 1 . 1 + 

h — — • m ~i 

p(3 1-/3/3 



/a 
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1 - /3 2 /? 2 



2.2 BFKL-like gluon polarization 



Now we take into account the BFKL-like gluon polarization 
the projector Pbfkl can be represented as 



As we shown in [E6 



>/3 

BFKL 



1 1 

2^ 



_ Ubx 



Q 2 



(13) 



In the previous subsection we have already presented the contributions to hard SF using 
the first term in the brackets of the r.h.s. of (0). Repeating the above calculations with the 
projector ~ q a q 13 , we obtain the total contributions to hard SF which can be represented as 
the following shift of the results in Eqs.([5[) - (|TT|): 



Ci(x) 


^2,BFKL { 


x), 


Cl(x) 


~ > Cl,bfkl( x ) 


/<i> _ 


_> f(l) 

~* J BFKL — 


1 




- 3bx 2 py 




/( 2) _ 


_, f(2) 

~* J BFKL — 


i 




- 3bx 2 f^ 





(14) 
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where 



l-x(l + b) 



x 



-2[x(l-x(l + b))(l + b-2a) + ap 2 ) ■ / x 



x(l -x(l + &))(!- 2a) -/ 2 



/ (2) = 4 • /3 (1 - (1 + 6)x) 2 



2 - (1 + 2te 2 ) ■ h - bx 2 ■ j 2 



(15) 



(16) 



Notice that the general formulae are needed only to evaluate the charm contribution to 
structure functions F 2 and F^, i.e. _F| and F[. To evaluate the corresponding light-quark 
contributions, i.e. F l L , we can use m 2 = limit of above formulae. 

2.3 The case m 2 = 

When m 2 = the hard SF C 9 k (x) are defined by f {2 \ f (1) and f {2) (as in Eqs.(|), © 



and (|14T) ) which can be represented as 



/ (2) 



2 - ( 1 - 2z(l + 6) + 2x 2 (l + 6) 2 J ■ L0) 

1 - 26x 2 • Hfl) 
1 - 26x 2 • Hj3) 



kc(l + 6)(l - (l + 6)x) 
(l + 6)(l-z(l + 6)) 



bx 



8bx 2 



f {2 \ 



/( 2 ) = A(l-x{l + b)f 



3 - (1 + 2bx 2 ) ■ L0) 



(17) 

(18) 
(19) 

(20) 



and, thus, (see Eq. (PI)) 



54 _ ,(1) 

P Jbfkl 



(-2) (1- a;(l + 6)) 



2 l-2x(l + 6) + 



1 - z(l + 6) - 4x 3 6(l + b) + 



z 2 (l - 6) 2 
1 - x(l + 6) 



l-x{l + b) 



(21) 



54 f (2) 
P ' J BFKL 



lx{l-x(l + b)) 



1 + 6 - 186x(l - x(l + 6)) 



+ 2te (3 - 4x(l + 6) + 6fo 2 (l - z(l + 6))) ■ L0) 



(22) 



where 



L( ^) = i.lni±4 
/? 1-/3 
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3 Relations between Fl, F2 and derivation of F2 in the 
case of collinear approximation 

Another information about the SF Fl can be obtained in the collinear approximation 
(i.e. when k\ = 0) in the following way. 

In the framework of perturbative QCD, there is the possibility to connecting Fl to F 2 
and its derivation dF 2 jd In Q 2 due the fact that at small x the DIS structure functions 
depend only on two independent functions: the gluon distribution and singlet quark one (the 
nonsinglet quark density is negligible at small x), which in turn can be expressed in terms 
of measurable SF F 2 and its derivation dF 2 j d\wQ 2 . 

In this way by analogy with the case of the gluon distribution function (see [0, |33[ and 



references therein), the behavior of F L (x, Q 2 ) has been studied in pT9[|-[pT|| , using the HERA 



data and the method [3B] f\ of replacement of the Mellin convolution by ordinary 

products. Thus, the small x behavior of the SF Fl(x,Q 2 ) can be extracted directly from 
the measured values of F 2 (x, Q 2 ) and its derivative without a cumbersome procedure (see 
These extracted values of Fl may be well considered as new small x 'experimental 
data 7 of Fl. The relations can be violated by nonperturbative corrections like higher twist 



ones (see [B8l B9]), which can be large exactly in the case of SF Fl [25, 40 



Because fc^-factorization approach is one of popular nonperturbative approaches used at 
small x, it is very useful to compare its predictions with the results of p^|-||21|| based on the 
relations between SF F L (x,Q 2 ), F 2 (x,Q 2 ) and dF 2 (x, Q 2 )/d\nQ 2 . It is the main purpose of 
the study. 

The /c^-factorization approach relates strongly to Regge-like behavior of parton distribu- 
tions. So, we restrict our investigations to SF and parton distributions at the following form 
(hereafter a = q,g): 

fa(x,Q 2 )~F 2 (x,Q 2 )~ x - 5 ^ (23) 

Note that really the slopes of the sea quark and gluon distributions: 5 q and 8 g , respec- 
tively, and the slope 5p2 of SF F 2 are little different. The slopes have a familiar property 



S q < d~F2 < o~ g (see Refs. [[|T|]-[[46] and references therein). We will neglect, however, this 
difference and use in our investigations the experimental values of S(Q 2 ) = 5f2{Q 2 ) extracted 
by HI Collaboration (see [f|l[] and references therein). We note that the Q 2 -dependence 
is in very good agreement with perturbative QCD at Q 2 > 2 GeV 2 (see ||48|| ). Moreover, the 
values of the slope S(Q 2 ) are in agreement with recent phenomeno logical studies (see, for 
example, [§]) incorporating the next-to-leading corrections f|9j (see also |]50|) in the frame- 
work of BFKL approach. 



Thus, assuming the Regge-like behavior ( P3"| ) for the gluon distribution and F 2 (x,Q 2 ) at 
x~ s ^> 1 and using the NLO approximation for collinear coefficient functions and anomalous 

9 The method is based on previous investigations pi [3^]]. 

10 Now the preliminary ZEUS data for the slope d\nF2/dln(l/x) are available as some points on Figs. 
8 and 9 in Ref. J47| . Moreover, the new preliminary HI points have been presented on the Workshop 
DIS2002 (see |l7]]). Both the new points are shown quite similar properties to compare with HI data pi] . 
Unfortunately, tables of the ZEUS data and the new HI data are unavailable yet and, so, the points cannot 
be used here. 
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dimensions of Wilson operators, the following results for Fl(x, Q 2 ) has been obtained in 



F L (x,Q 2 



Bf +5 (l + a s (Q 2 )R y L 



+ 



(0),l+<5 
799 

a s (Q 2 ) (B 



a s (Q 2 



B 



i _(l),l+<5 

+ Tqg 

q,l+S 

g,i+8 >qg 



7, 



(0),l+<5 



'/'/ 



dF 2 (x£,Q 2 ) 
d\nQ 2 

F 2 (x£,Q 2 ) 



+ 0(a 2 s ,x 2 ' 



-^ax 1 - 8 ) 



(24) 



where 



<qg 



7$}'" + B 



!qq J 



R 



Rf 1 - B 9 2 ' n 



B 



q,v 



Here 7«>" and T?£ n 
of Wilson operators 7?? 



= q, g) are the combinations of the 'anomalous dimensions' 
0,7$)'" + a^W'" + O(a^) and 'Wilson coefficients' a.fl^l + 

a^i?^) + O(a^) and a s B 2 ,v + 0(a 2 ) with the 'moment' argument 77 (i.e., the combinations 
of the functions which can be obtained by analytical continuation of the corresponding 
anomalous dimensions and coefficient functions from integer values n of their argument to 
non- integer ones n). 

Note that, in principle, any term like ~ l/(n + m) (m = 0, 1,2, ...) which comes to the 
corresponding combinations of the anomalous dimensions and coefficient functions: 7^' ra 
and R a f n , should contribute to Eq.(p4[) in the following form (after replacement of Mellin 
convolutions by usual products in the DGLAP equations (see p6|)): 



1 + 5 + m 



1 + 



T(2 + 5 + m)T(l + u) 
T(2 + 5 + m + 1 + v) 



x 



1+5H 



(25) 



where the value of v comes (see [||T], [52|, [53|]) from asymptotics of parton distributions f a (x) 



at x 
term 



1: fa 



x 



X 



1+d + in 



in Eq. 



and v » 4 from quark account rules |54|. The additional 



is important only for rn — —1 case (i.e. for the singular parts 



~ l/(n — 1) of the corresponding anomalous dimensions and coefficient functions) and quite 
small values of 5 (i.e. for x 5 ~ Const). 

Thus, excepting the case when m = — 1 and x 5 ~ Const, we can replace (^) by its 
first term 1/(1 + 5 + m), i.e. our variables Jy' 11 and R£ are just the combinations of 
corresponding anomalous dimensions and coefficient functions at n — 1 + 8. When m — — 1 
and x s ~ Const at the considered small x range, we should replace the term l/{n — 1), if it 
was exist in the variables 7^ and Rf (a = g, g), by the following term: 



/ qa 
1 



r(i - 6)v(i + u) j 
m-6 + u) 



(26) 



Note also that the 1/5 coincides approximately with 1/5 when 5^0 and x —>■ 0. However, 
at 5 — > 0, the value of 1/5 is not singular: 



i In fi 

5 \x 



(27) 



11 Bccause we consider here Fa (a;, Q 2 ) but not the singlet quark distribution in the corresponding DGLAP 
equations ]3C|]. 



12T 



In our formula (24) we have our interest mostly to gluons, so we can apply v = v g w 4 below. 
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So, the Eq. 
a, b 



and 



4|) together with well-known expressions of anomalous dimensions 7„ b 

a = q,g) 



(0),n 



g, g) and coefficient functions B 2 and 



sec 



55! and 



56| , |57| , respectively, and references therein) gives a possibility to extract SF Fl at small x 



values. The calculations are based on precise experimental data of SF F 2 and its derivatives 
dF 2 /d\n(Q 2 ) and 5 = d\nF 2 /d\n(l/x). 

For concrete 5- values, the Eq. fl24j ) simplifies essentially (see [p0|| ). For example, for 
0.3 we obtain (for the number of active quarks / = 4 and MS 1 scheme): 



Fl(z,Q 



0.84 



+ 



l + 59.3a,(Q 2 ) 



rfF 2 (0.48x,Q 2 
dlnQ 2 



+ 3.59a s (Q 2 )F 2 (0.48x,Q 2 



(28) 



At arbitrary 5 values, in real applications it is very useful to simplify Eq. ( |2"4| ) as follows. 
We keep the exact ^-dependence only for the leading order terms, which are very simple. In 
the NLO corrections we extract the terms ~ 1/5, which are changed strongly when < 5 < 1, 
and parameterize the rest terms in the form: + b^5 + c-fi 2 . The coefficients a^, 6j, Cj are 
fixed from the agreement these parameterizations with the exact values of 7^ ,?? and R L ' at 
5 = 0, 0.3 and 0.5. These exact values can be found in Refs. [[33|, [2(], |2T|j . 

Then, the approximate representation of Eq. fl2"4|) of for arbitrary 5 value has the form: 



F L {x,Q 2 



r(i+5)mr 



'l + 30a s (Q 2 )[l/5- 1 -§p 1 (5)}) 



dF 2 (x£(5),Q 2 ) 8 



dlnQ 2 



+ - P 2(5)a s (Q 2 )F 2 (x£(5),Q 2 



+ 0(a 2 ,a s x,x z 



where 



(29) 



r{5) 
Pi{5) 



15 



2 + 5 + 5 2 ' 
1 + 5 + 5 2 /A. 



m = 

P2(5) 



r(5) 



r(l + S)' 

= 1 - 2.395 + 2.69<5 2 , 



(30) 



4 Comparison with Fl experimental data 



With the help of the results obtained in the previous Section we have analyzed experi- 
mental data for SF F L Q from HI |T3j and |§, NMC @, CCFR [§ and [§, BCDMS 



6TJ1 collaborations 0. Note, that we do not correct CCFR data [59 and [^0| which has been 
obtained in uN processes because the terms ~ x ■ m 2 /Q 2 , which are different in pN and uN 
processes, are not so strong at low x. 

We calculate the SF Fl as the sum of two types of contributions: charm quark one F£ 
and light quark one F l L : 



Ft, + Fl 



13 Sometimes there are experimental data for the ration R = 
F L because F L = F 2 R/(l + R). 

14 We do not use experimental data of R from SLAC pi 
because they are obtained at quite large x values. 



(31) 

cl/ct, which can be recalculated for the SF 
1, EM ii and CDHSW 161 Collaborations 
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We use the expression (|3|) for the calculation of the both SF F L and F£ in the following form 
(here a c = m 2 c /Q 2 ): 



F l L (x,Q 2 ) 



E 



f=uds 

2 r z (l) 

1 = 1 mm 



ld ^Cl iBFKL (- Q 2 1 0)zf g (z,Ql) + 



z 



, X 



2(0 



dk ± Cf BFKL (—, Q , k ± ) $(z, k ± , Q ) 



F c L (x,Q 2 ) 



) -Cl,BFKL(-,mlQ 2 ,0)zf g (z J Ql) + 

x(l+4a c ) Z Z 



+ E 



,(«) J „L.2(i) 



(0 



dfci Cl BFKL (- z , ml, Q 2 , k\) 3>(z, kl, Q, 



where C 9 L bfkl( x b, Q 2 , m>l, kj_) are given by Eqs. (^) and (|T^). 

The integration limits in the expression (K3|) have the following values: 



z {1) 

"■'mm 

,2(1) 
h -Lmin 

J 2 ) 



V ' i Q 2 ' maa 

O 2 k 2(l) 



2x(l + 2a c ); 



A: 



2(2) 
_Lmin 



2a;(l + 2a c ), z®, 
O 2 k 2{2) 



(--(l + 4a c ))Q 



2. 



1: 



Q 2 ; 



(32) 



(33) 



(34) 



The ranges of integration correspond to positive values of square roots in expressions 
(PH), QTTD, (|15|) and (|16|) and also should obey to kinematical restriction (z < (l + 4a c + fe) -1 ) 
following from condition [3 2 > (see Eqs. ([I0|)-(|l2"l)). In Eq.(^) the ranges (|34]) are used at 
a, = 0. 



In Fig. 2 we show the SF Fl as a function x for different values of Q 2 in comparison 



with HI experimental data sets: old one of |13| (black triangles), last year one of |15 
(black squares) and new preliminary one of jn| (black circles) and also with NMC 



(white triangles), CCFR |59|] (white circles) and BCDMS |6l] data (white squares). For 
comparison with these data we present the results of the calculation with three different 
parameterizations for the unintegrated gluon distribution $(x, kj_, Qq) at Qq = 4 GeV 2 . All 



of them: Kwiecinski-Martin-Stasto (KMS) one JT2] and Blumlein (JB) one |68 and Golec- 
Biernat and Wusthoff (GBW), have been used already in our previous work PB1 and reviewed 
there. 

There are several other popular parameterizations (see, for example, Kimber-Martin- 
Ryskin (KMR) J70[ and Jung-Salam (JS) f7lf), which are not used in our study mostly 
because of technical difficulties. Note that all above parameterizations give quite similar 
results excepting, perhaps, the contributions from the small fc^-range: k\ < 1 GeV 2 (see 



Ref. (|TTJ and references therein). Because we use Ql = 4 GeV 2 in the study of SF F L , our 
results depend very slightly on the the small /c^-range of the parameterizations. In the case 
JB, GBW and KMS sets this observation is supported below by our results and we expect 
that the application of KMR and JS sets should not strongly change our results. 

The differences observed between the curves 2, 3 and 4 are due to the different behavior 
of the unintegrated gluon distribution as function x and k±. We see that the SF F L obtained 
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in the /^-factorization approach with KMS and JB parameterizations is close each other Q 
and higher than the SF obtained in the pure perturbative QCD with the GRV gluon density 
at the leading order approximation. Otherwise, the /^-factorization approach with GBW 
parameterization is very close to pure QCD predictions Q: it should be so because GBW 
model has deviations from perturbative QCD only at quite low Q 2 values. Thus, the pre- 
dictions of perturbative QCD and ones based on kx factorization approach are in agreement 
each other and with all data within modern experimental uncertainties. So, a possible high 
values of high-twist corrections to SF Fl predicted in can be important only at low Q 2 
values: Q 2 < Q 2 = 4 GeV 2 . 

Fig. 3 is similar to Fig. 2 with one excepting: we add 'experimental data' obtained using 
the relation between SF Fl(x, Q 2 ), F 2 (x, Q 2 ) and dF 2 (x, Q 2 ) / 'dliiQ 2 (see Section 3) as black 
stars. Because the corresponding data for SF F 2 (x, Q 2 ) and dF 2 (x, Q 2 )/d In Q 2 essentially 
more precise (see [|l5j) to compare with the preliminary data |18| for Fl, the 'experimental 
data' have strongly suppressed uncertainties. As it is shown on Fig. 3 there are very good 
agreement between the new preliminary data fl~8| , the 'experimental data' and predictions 



of perturbative QCD and /c^-factorization approach. 



To estimate of the value of charm mass effect, we recalculate SF F£ also in massless 
approximation similar to fl3"2"|). In Fig. 4 we show importance of exact m c -dependence in 
hard SF of F[ to compare with its massless approximation, where we should have Fl(m c = 
0)/Fi(m c = 0) = 2/5. As it is possible to see on Fig. 4, the ratio F[/Fl goes to massless 
limit 2/5 only at asymptotically large Q 2 values. 



5 Conclusions 

We have applied in the framework of /^-factorization approach the results of the cal- 
culation of the perturbative parts for the structure functions Fl and F[ for a gluon target 
having nonzero momentum square, in the process of photon-gluon fusion to the analysis of 
present data for the structure function Fl [_} The analysis has been performed with several 
parameterizations of unintegrated gluon distributions, for comparison. We have found good 
agreement between all existing experimental data, the predictions for Fl obtained from the 
relation between SF Fl(x, Q 2 ), F 2 (x, Q 2 ) and dF 2 (x, Q 2 )/d\nQ 2 and the results obtained in 
the framework of perturbative QCD and ones based on /^-factorization approach with the 
three different parameterizations of unintegrated gluon distributions. 

We note that it could be also very useful to evaluate the SF F 2 itself and the derivatives 
of F 2 respect to the logarithms of 1/x and Q 2 with our expressions using the unintegtated 
gluons. We are considering to present this work and also the predictions for the ratio 
R = ctl/(Tt in a forthcoming article. 

15 Note that very similar results have been obtained also for Ryskin-Shabelsky parameterization |72| (see 
f§). 

16 This fact is evident also from quite large value of Q\ — 4 GeV 2 chosen here. 

17 In Ref [^6| we have also obtained quite large contribution of SF F£ at low x and high Q 2 (Q 2 > 30 
GeV 2 ). 

18 A study of the SF F2 in the framework of fcT-factorization has been already done in Ref. J74j . 
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The consideration of the SF F 2 in the framework of the leading-twist approximation of 
perturbative QCD (i.e. for "pure" perturbative QCD) leads to very good agreement (see 
Ref. f|5|| and references therein) with the HERA data at low x and Q 2 > 1.5 GeV 2 . The 



agreement improves at lower Q 2 when higher twist terms are taken into account in Ref. 
| 38| , |39[| . As it has been studied in Refs. j|5], |38j , the SF F 2 at low Q 2 is sensentive to the 



small-x behavior of quark distributions. Thus, our future analysis of F 2 in broad Q 2 range in 
the framework of /^-factorization approach should require the incorporation of parameteri- 
zations of unintegrated quark densities, introduced recently (see Ref. |7[], |TTJ and references 
therein). 
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Figure captions 



Fig. 1 The diagrams contributing to T^ v for a gluon target. They should be multiplied by 
a factor of 2 because of the opposite direction of the fermion loop. The diagram (a) should 
be also doubled because of crossing symmetry. 



Fig. 2 The structure function Fl(x,Q 2 ) as a function of x for different values of Q 2 com- 
pared to experimental data. The HI data: the first 1997 ones [13|], new 2001 ones |H| and 
preliminary ones |18| are shown as black triangles, circles and squares, respectively. The 
data of NM |55 |, CCFR 0] and BCDMS Collaborations are shown as white triangles, 
circles and squares, respectively. Curves 1, 2, 3 and 4 correspond to SF obtained in the 
perturbative QCD with the GRV [^31 gluon density at the leading order approximation and 
to SF obtained in the factorization approach with JB (at Q\ = 4 GeV 2 ) [68|, Kwiecinski- 



Martin-Stasto (KMS) and GBW 591 parametrizations of unintegrated gluon distribution. 



Fig. 3 The structure function Fl(x,Q 2 ) as a function of x for different values of Q 2 . To 
compare with Fig. 2 'experimental data' (see fll9"l-pl] and Section 3) are added as black 
stars. The 'experimental data' values depend mostly on the derivative dF 2 (x, Q 2 )/d In Q 2 , 
which data are know at little different Q 2 values (see ]TB|]). So, the 'experimental data' 
obtained at 12, 15, 20, 25 and 35 GeV 2 are presented here at 13.4, 15.3, 22.4, 29.6 and 39.7 
GeV 2 , respectively. 

Fig. 4 The ratio F£(x, Q 2 )/Fl(x, Q 2 ) as a function of x for different values of Q 2 . 
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Figure 1: 
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